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Abstract. I discuss recent results on lattice QCD calculations with the main emphasis on the 
thermodynamics of the crossover region, restoration of the chiral symmetry and fluctuations 
of conserved charges as indicator of deconfinement, that may also be used to determine the 
chemical freeze-out conditions in heavy-ion collision experiments. 



1. Introduction 

It has been established by lattice QCD calculations that at the physical values of light quark 
masses and at vanishing chemical potential there is no genuine phase transition in QCD, but 
rather a "rapid" crossover [H El E] • While the low-temperature phase exhibits confinement and 
breaking of chiral symmetry, at high temperatures the behavior of the theory is qualitatively 
different - the interaction between quarks and gluons decreases due to asymptotic freedom, 
leading to deconfinement, and the chiral symmetry is restored (see [I] for a recent review). In 
the crossover region the QCD partition function does not exhibit a singularity, so it is not a 
surprise that different physical observables show a change in their temperature dependence at 
somewhat different temperatures. 

2. Restoration of chiral symmetry and T pc in QCD 

Let us start the discussion with the phenomenon of chiral symmetry restoration. The QCD 
Lagrangian at zero light quark mass allows for independent left and right rotations, i.e. possesses 
a SU(2) l x SU (2)r symmetry. This symmetry is however spontaneously broken in the vacuum, 
as indicated by a condensate of quark anti-quark pairs. On general grounds [5J it is expected 
that at zero light quark mass there is a second-order phase transition in QCD (in the 0(4) 
universality class) for which the chiral condensate is the order parameteiQ. The renormalized 
quark condensate A/ s (I reserve the name "chiral condensate" for the light quark condensate at 
zero mass) is shown in Fig. [1] and defined as: 

\s = -rrr\ ar ; i , ; » = tt^— = ttttt (™„ ), q = l,s, (1) 
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This is true under an additional assumption that the U(1)a axial symmetry breaking is still large at the chiral 
critical temperature, T c , and the symmetry is restored at higher temperatures. Otherwise the symmetry group 
would be enlarged, and the 0(4) universality class would be no longer appropriate. 



where M q is the fermion matrix, i.e. the discretized version of the Dirac operator, N^N T is 
the four-dimensional lattice volume. The disconnected susceptibility which is proportional to 
fluctuations of the quark condensate 

XP C = {((TrMrf) - (TrM-) 2 } (2) 

is shown in Fig. [2j The results are obtained with the asqtad and HISQ/tree action at the light to 
strange quark mass ratio mi = m s /20 (somewhat heavier but close to the physical mi = m s /27). 
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Figure 1. The renormalized quark condensate 
with the asqtad and HISQ/tree action. 
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Figure 2. The disconnected susceptibility 
with the asqtad and HISQ/tree action. 



If the physical light quark mass is small enough one may hope that remnants of criticality (plus 
subleading corrections) still govern the crossover region. Applicability of the critical scaling was 
extensively studied in Ref. [6j with staggered fermions and the scaling was indeed observed. (An 
additional complication for staggered fermions is that at non-zero lattice spacing the relevant 
universality class is 0(2) rather than 0(4). Luckily, the numerical difference between 0(2) 
and 0(4) fits is not significant for our discussion, comparable with other uncertainties, and 
we refer to this analysis as O(N) scaling.) This approach has been taken by the HotQCD 
collaboration, who calculated the light quark condensate and its susceptibility with the asqtad 
and HISQ/tree actions on lattices with the temporal extent N T = 6, 8 and 12 at several values 
of the light quark mass and then performed fits to O(N) scaling functions complemented by 
non-singular terms. The pseudo-critical temperature, T pc , defined this way as a location of the 
peak of the chiral susceptibility, reduces to the true critical temperature in the chiral limit. The 
result in the continuum limit at the physical light quark mass is T pc = 154(9) MeV [7j. This 
is compatible with earlier results by the Budapest- Wuppertal collaboration that are in a range 
T pc = 147 — 157 MeV, depending on what chiral observable is picked to determine T pc [8]. 

The studies discussed so far rely on staggered fermions (one of the numerically cheapest lattice 
fermion formulations). The fermion doubling problem leads to 16 species of lattice fermions per 
one in continuum and staggered fermions reduce them to 4 species. The latter are however 
non-degenerate (and are labeled by "taste" to distinguish from flavor). In other words, taste 
symmetry is broken, for particular actions we use at 0(a 2 ), a being the lattice spacing. This leads 
to a distorted hadron spectrum, and taste-breaking effects have been identified as the largest 
source of systematic errors in staggered simulations (a discussion in relation to thermodynamics 
is presented in Ref. [9]). It is desirable to have independent tests of the results also with 
other types of lattice fermions. The latter require more computational resources. Calculations 
directly at the physical mass thus may not be feasible. However, some conclusions can be drawn 



from simulations at heavier masses. The Budapest- Wuppertal collaboration showed agreement 
between the continuum-extrapolated light quark condensate for staggered and Wilson fermions 
at the pion mass = 545 MeV [10J, Fig. [3] (left). Another study by the same group compared 
the condensate with overlap fermions (that provide exact chiral symmetry at non-zero lattice 
spacing) again to staggered calculations at m w = 350 MeV [11]. Although the continuum 
limit for overlap fermions has not been taken, it seems, Fig. [3] (middle), that disagreement 
would be unlikely. The HotQCD collaboration pursued simulations with domain-wall fermions 
at m n = 200 MeV [T2]. A comparison of the staggered and domain- wall disconnected chiral 
susceptibility is shown in Fig. [3J (right). Note similar location of the peak in the domain- wall 
(DWF) and HISQ/tree N T = 12 data. (Strictly speaking, results for different lattice actions 
should be compared in the continuum limit, however, Fig. [3] (right) provides a very encouraging 
check.) 
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Figure 3. Comparison of the observables indicating restoration of the chiral symmetry between 
staggered and three other fermion formulations (Wilson, overlap and domain- wall) . See text for 
details. 



3. Deconfinement and fluctuations of conserved charges 

QCD with infinitely heavy quarks, i.e., SU(3) pure gauge theory, exhibits a first-order phase 
transition at T c ~ 260 MeV [13] . The order parameter is the renormalized Polyakov loop: 

L ren (T) = exp(-cN T /2) ^Tr j| U (x , x)^ , (3) 

where Uq(xo,x) is the gauge field in the time direction. The exponential prefactor takes care of 
the multiplicative renormalization. L ren can be related to the free energy Fq(T) of a static test 
quark as L ren (T) = exp(— Fq/T). In pure gauge theory the Svetitsky-Yaffe argument [H] allows 
to interpret deconfinement as a transition from the symmetric phase, L ren = 0, into the broken 
phase, L ren ^ 0, where the Polyakov loop assumes one of the three equally possible values, 
related by Z3 symmetry. So, historically, based on the experience from pure gauge theory, the 
Polyakov loop served as an indicator of deconfinement in lattice QCD. Dialing quark masses to 
smaller values breaks the Z3 symmetry (making one of the values of L ren preferable over others, 
just like a magnetic field in spin systems). Due to string breaking by light dynamical quarks in 
QCD the value of L ren is always non-zero in the confined phase. Thus, in QCD the Polyakov 
loop loses its meaning of the order parameter and is not related to any singularity of the partition 
function (in the chiral limit). As one can see in Fig. 0] (left) the temperature dependence of the 
Polyakov loop in SU(2) and SU(3) pure gauge theory and in QCD is quite different. Sharp 
behavior in the former case gives way to smoother change over wide temperature range. The 
behavior of the Polyakov loop in QCD at low temperatures can be understood in terms of the 
hadronic degrees of freedom |15|. I16|. However, this description breaks down in the vicinity of 
the transition region. Its role as indicator for deconfinement thus is obscured. 

Another possibility to probe deconfinement is through fluctuations of conserved charges. For 
instance, light and strange quark number fluctuations, that are suppressed at low temperature, 
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Figure 4. Temperature dependence of the Polyakov loop in SU(2) and SU (3) pure gauge 
theory and QCD (left). Comparison of the behavior of the renormalized light quark condensate 
and the Polyakov loop (middle) and light and strange quark number fluctuations (right). 

signal liberation of the degrees of freedom with quantum numbers of quarks. In Fig. 0] (middle, 
right) we compare the light quark condensate to the Polyakov loop, light and strange quark 
number susceptibility. As one can see, the observables associated with deconfinement rise 
somewhat slowly and over substantial temperature range. This reflects that singular terms 
contributing to quadratic fluctuations are still subleading and only the dominant temperature 
dependence of regular terms is seen. The analysis of the temperature dependence of higher order, 
e.g. fourth order, susceptibilities is needed to become sensitive to deconfining features. Thus 
in full QCD associating a particular temperature with deconfinement by searching for inflection 
points or alike in the Polyakov loop or quadratic fluctuations does not appear meaningful. 
Rather, full temperature dependence of various observables should be studied. In particular, one 
can study fluctuations and correlations of baryon number B, electric charge Q and strangeness 
S that are given by the derivatives of the pressure with respect to the corresponding chemical 
potentials 

BQS = g^g) 

X » k d(fiB /TY d{^Q /T)i d{ns /T) k ' [ 1 

The fluctuations of strangeness, baryon number and electric charge, calculated with the 
HISQ/tree action and extrapolated to the continuum, are shown in Fig. [5] (HotQCD 
collaboration [17J). Solid curves represent the Hadron Resonance Gas model, which is in 
remarkable agreement with the lattice data up to T ~ 150 — 160 MeV. (Electric charge 
fluctuations are the most sensitive to the cutoff effects in the pion sector and obtaining the 
continuum limit is, thus, a more subtle issue.) The results from the Budapest- Wuppertal 
collaboration |18| for the same quantities are in complete agreement. 
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Figure 5. Fluctuations of strangeness, baryon number and net electric charge with the 
HISQ/tree action extrapolated to the continuum and compared to the Hadron Resonance Gas 
model. 



To probe the critical behavior of QCD further one can consider higher-order cumulants of 
conserved charges. Comparing ratios of cumulants calculated on the lattice with experimental 



data, under assumption that at the time of chemical freeze-out the system created in heavy-ion 
collisions can be described by equilibrium thermodynamics with temperature Tf and baryon 
chemical potential \J B , one can extract these freeze-out parameters. A recent attempt to use 
electric charge fluctuations for this purpose has been presented in Ref. [19]. By constraining the 
mean net strangeness and electric charge to those of the incident nuclei, 



M s = 0, M Q = rM B 



0.4 in gold-gold and lead-lead collisions), 



(5) 



one can determine the strange and electric charge chemical potential of the system, 
next-to-leading order (NLO) in the baryon chemical potential: 
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The qi, s, coefficients are related to ratios of generalized susceptibilities, i.e. various combinations 
of derivatives of the partition function with respect to chemical potentials, Eq. @. It has been 
shown in Ref. [TH] that NLO corrections are small and at NLO the dependence of [IqI\xb and 
Ms/Mb 011 Ms at several temperatures is shown in Fig. [UJ The bands give total uncertainty of 
the lattice calculation and dashed lines represent the HRG model. Once fiQ and fig satisfying 
the constraint (0) are fixed, ratios of cumulants can be evaluated, and Ref. [TO] focused on 
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(8) 



where X = B or Q, Mx is mean, a\ is variance and Sx skewness for corresponding conserved 
charges. The results are shown in Fig.[7|and[8l If is determined from experiment, one can use 
the band in Fig. [7] to convert it to the freeze-out temperature Tf. Then from an experimentally 
determined ratio Rf 2 , Tf and the corresponding band in Fig. [8] one can find the freeze-out 
chemical potential fj, B . Once the full experimental analysis of R^ and R® 2 is available, the 
results of Ref. |19j can determine Tf and \x B for various center-of-mass energies and thus map a 
part of the freeze-out curve on the QCD phase diagram. 
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4. A list of other results on the lattice 

Due to the space constraints many other recent interesting calculations on the lattice could not 
be included. Among those I would like to mention: 

• Progress on the equation of state with 2+1 flavors of staggered fermions: an update of the 
previous results from the Budapest- Wuppertal [20] and HotQCD [21] [22] Collaboration. 
Preliminary results still indicate a discrepancy in 150 — 300 MeV range, which may be 
resolved by ongoing calculations. 

• 2+1 flavor equation of state with Wilson fermions from the WHOT-QCD collaboration |23j. 

• The equation of state with dynamical charm quark, see Refs. [20J and [24j for preliminary 
results. 

• Calculation of the charmonium screening masses from spatial correlation functions in 2+1 
flavor QCD [25] and Maximum Entropy Method (MEM) reconstruction of the charmonium 
spectral functions in quenched QCD |26j. 

• Thermodynamics of the crossover region with high magnetic fields [27]. 
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